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Four-dimensional cylindrically symmetric spacetimes with homothetic self-similarity are studied 
in the context of Einstein's Theory of Gravity, and a class of exact solutions to the Einstein-massless 
scalar field equations is found. Their local and global properties are investigated and found that 
they represent gravitational collapse of a massless scalar field. In some cases the collapse forms 
black holes with cylindrical symmetry, while in the other cases it does not. The linear perturbations 
of these solutions are also studied and given in closed form. From the spectra of the unstable eigen- 
\ modes, it is found that there exists one solution that has precisely one unstable mode, which may 

. represent a critical solution, sitting on a boundary that separates two different basins of attraction 

' in the phase space. 
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I. INTRODUCTION 



m 

^ ' The studies of non-linearity of the Einstein field equations near the threshold of black hole formation reveal very 
rich phenomena [1] , which are quite similar to critical phenomena in Statistical Mechanics and Quantum Field Theory 
I [2]. In particular, by numerically studying the gravitational collapse of a massless scalar field in 3 -f 1-dimensional 
spherically symmetric spacetimes, Choptuik found that the mass of such formed black holes takes the form. 



Mbh^C{p){p~p*)\ (1.1) 

where C (p) is a constant and depends on the initial data, and p parameterizes a family of initial data in such a way 
that when p > p* black holes are formed, and when p < p* no black holes are formed. It was shown that, in contrast 



O 

I , to C{p), the exponent 7 is universal to all the families of initial data studied, and was numerically determined as 
. ry ^ 0.37. The solution with p = p* , usually called the critical solution, is found also universal. Moreover, for the 
rrl massless scalar field it is periodic, too. Universality of the critical solution and the exponent 7, as well as the power-law 
^ • scaling of the black hole mass all have given rise to the name Critical Phenomena in Gravitational Collapse. 
k>( \ Choptuik's studies were soon generalized to other matter fields [3]. From all the work done so far, the following 
. seems clear: (a) There are two types of critical collapse, depending on whether the black hole mass takes the scaling 
' form (1.1) or not. When it takes the form, the corresponding collapse is called Type // collapse, and when it does 
not it is called Type / collapse. In the type // collapse, all the critical solutions found so far have either discrete 
self-similarity (DSS) or homothetic self-similarity (HSS), depending on the matter fields. In the type / collapse, the 
critical solutions have neither DSS nor HSS. For certain matter fields, these two types of collapse can co-exist, (b) For 
Type // collapse, the corresponding exponent is universal only with respect to certain matter fields. Usually, different 
matter fields have different critical solutions and different exponents. But for a given matter field the critical solution 
and the exponent are universal ^ . (c) A critical solution for both of the two types has one and only one unstable mode. 
This now is considered as one of the main criteria for a solution to be critical, (d) The universality of the exponent is 
closely related to the number of unstable modes. In fact, the unstable mode, say, fci, of the critical solution is related 
to the exponent 7 via the relation, 7 = |fci| ^, which can be obtained by using dimensional analysis [4]. 



*E-mail: Anzhong_Wang@baylor.edu. On leave from Department of Theoretical Physics, the State University of Rio de 
Janeiro, Brazil. 

^So far, the studies have been mainly restricted to spherically symmetric case and their non-spherical linear perturbations 
[3]. Therefore, it is not really clear whether or not the critical solution and exponent are universal with respect to different 
symmetries of the spacetimes. 
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From the above, one can see that to study critical collapse, one may first find some particular solutions by imposing 
certain symmetries, such as, DSS or HSS. This can simplify the problem considerably. For example, in the spherically 
symmetric case, by imposing HSS symmetry the Einstein field equations will be reduced from PDE's to ODE's. 
Once the particular solutions are known, one can study their linear perturbations and find out the spectrum of the 
corresponding eigen-modes. If a solution has precisely one unstable mode, it may represent a critical solution, sitting 
on a boundary that separates two different basins of attraction in the phase space. 

The studies of critical collapse have been mainly numerical so far, and analytical ones are still highly hindered by 
the complexity of the problem, even after imposing some symmetries. Lately, some progress has been achieved in the 
studies of critical collapse of a masslcss scalar field in an anti-de Sitter background in 2 + 1-dimensional spacetimes 
both numerically [5,6] and analytically [7-9]. This serves as the first analytical model in critical collapse. 

In this paper, we shall present another analytical model that represents critical collapse of a massless scalar field 
in four-dimensional Einstein's Theory of Gravity with cylindrical symmetry. Although spacetimes with cylindrical 
symmetry do not represent realistical models, the studies of them can provide deep insight into the nonlinearity of 
the Einstein field equations. In particular, they may shine some light on the possible roles that gravitational radiation 
and angular momentum may play in critical collapse. In fact, such studies have already shown to be very useful in 
probing non-spherical gravitational collapse [10]. In addition, they may also provide a useful testbed for numerical 
relativity [11] and Quantum Gravity [12]. 

The rest of the paper is organized as follows: In Sec. II we first review the regularity conditions for a four-dimensional 
cylindrical spacetime, including the ones at the symmetry axis. Then we introduce the notion of homothetic self- 
similarity with cylindrical symmetry. In Sec. Ill, a class of exact solutions with such a symmetry to the Einstcin- 
massless scalar field equations is presented. It is shown that they represent gravitational collapse of a scalar field, 
in which black holes can be formed. In Sec. IV, the linear perturbations of these solutions are studied and given in 
closed form. After properly imposing boundary conditions, the spectra of the unstable modes of the perturbations are 
determined. In particular, it is found that there exists a solution that has precisely one unstable mode, which may 
represent a critical solution, sitting on a boundary that separates two different basins of attraction in the phase space. 
In Sec. V, the main results are summarized and some concluding remarks are given. There are also two appendices, 
A and B. In Appendix A, the Ricci tensor is given in terms of self-similar variables. The linear terms of perturbations 
of the Ricci tensor are also given there. In Appendix B, the expansions of the out- and in-going radial null geodesies 
are calculated, from which trapped surfaces and apparent horizons are defined. 



II. SPACETIMES WITH HOMOTHETIC SELF-SIMILARITY 



The general metric for cylindrical spacetimes with two hypersurface orthogonal Killing vectors takes the form [13], 

ds^ = e-^^*'--) [dt^ - dr^) - r^e-^f*''-) (e^^'^^Uw^ + e-^^''^Ue^) , (2.1) 

where = {t, r, w, 9} are the usual cylindrical coordinates, and the hypersurfaces 6 = 0,2tt are identified. The two 
Killing vectors are given by ^(^y) = and ^(g) = dg. To have cylindrical symmetry, some physical and geometrical 
conditions needed to be imposed. In general this is not trivial. As a matter of fact, when the symmetry axis is 
singular, it is still an open question: which conditions should be imposed [14]. Since in this paper we are mainly 
interested in gravitational collapse, we would like to have the axis regular in the beginning of the collapse. By this 
way, we are sure that the singularity to be formed later on the axis is indeed due to the collapse. Thus, following [15] 
we impose the following conditions: 

(i) There must exist a symmetry axis. This can be written as 



X 



^(6)^(0)9^^ ^0, (2.2) 



as r — > 0"*", where we have chosen the radial coordinate r such that the axis is located at r = 0. 
(ii) The spacetime near the symmetry axis is locally flat. This can be expressed as [13] 

as r ^ 0+, where ( )_o, = 9( )/dx". Note that solutions failing to satisfy this condition are sometimes acceptable, 
and are usually expected that the singularities located on the axis should be replaced by some kind of sources in 
more realistic models. A particular case of these is when the right-hand side of the above equation approaches a finite 
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constant, and the singularity now can be related to a line-like source [16]. In this paper, since we are mainly interested 
in gravitational collapse, we shall not consider these possibilities and assume that the above condition holds strictly 
at the initial of the collapse. 

(iii) No closed timelike curves (CTC's). In spacetimes with cylindrical symmetry, CTC's can be easily introduced. 
To guarantee their absence, we impose the condition 

^(0)^10)9,.^ < 0, (2.4) 

in the whole spacetime. 

In addition to these conditions, it is usually also required that the spacetime be asymptotically flat in the radial 
direction. However, since we consider solutions with self-similarity, this condition cannot be satisfied by such solutions, 
unless we restrict the validity of them only up to a maximal radius, say, r = ro{t), and then join the solutions with 
others in the region r > ro{t), which are asymptotically flat in the radial direction. In this paper, we shall not consider 
such a possibility, and simply assume that the self-similar solutions are valid in the whole spacetime. 

Spacetimes with homothetic self- similarity (or self- similarity of the first kind) is usually defined by the existence of 
a conform Killing vector ^'^ that satisfies the equations [17], 

where a semicolon ";" denotes the covariant derivative. It can be shown that for the spacetimes given by Eq.(2.1) the 
conditions (2.5) imply that 

M{t,r)=M{z), S{t,r) = S{z), V{t,r) = V{z), (2.6) 
where the self-similar variable z and the corresponding conform Killing vector ^'^ are given by 

ed^=tdt + rdr, ^=^- (2.7) 

It is interesting to note that under the coordinate transformations 

t = ait-\-a2f, r = ast-\-a4f, (2.8) 

the metric (2.1), the regular conditions (2.2)-(2.4), and the self-similar conditions (2.6) and (2.7) are all invariant, 
where a,'s are real constants, subject to aia2 — 0304 = 0. Using this gauge freedom, we shall assume that 

M{t, 0) = 0, (2.9) 

that is, the timelike coordinate t measures the proper time on the axis. 



III. SELF-SIMILAR SOLUTIONS OF MASSLESS SCALAR FIELD 



For a massless scalar field, the Einstein field equations read 

-R/ni/ = i^4>,ti4>,u, (3.1) 

where n[= SttG/c^] is the Einstein coupling constant. In this paper we shall choose units such that k = 1. The scalar 
field satisfies the Klein-Gordon equation, 

g''^c^:»l3 = 0. (3.2) 

However, this equation is not independent of the Einstein field equations (3.1) and can be obtained from the Biachi 
identities G^a-f^g"^ = 0. 

On the other hand, it can be shown that a massless scalar field 4>{t,r) that is consistent with spacetimes with 
homothetic self-similarity must take the form, 

ct>{t,r)=2qln{-t)-i-<fi{z), (3.3) 

where q is an arbitrary constant, and ip{z) is a function of z only, which will be determined by the Einstein field 
equations (3.1). Inserting Eqs.(A.2) and (3.3) into Eq.(3.1) and considering the self-similar conditions (2.6), we find 
the following solutions 
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M{z) = 2q^ In (l - z"^) , S{z) = ln{z), 

Viz) = -ln{z), ip{z) = 0. (3.4) 

When q = the corresponding spacctime is flat. Thiis, in the following we shall assume that 17 7^ 0. Then, it 
can be shown that these solutions satisfy all the conditions (2.2)-(2.4) and (2.9), and the corresponding Ricci and 
Kretschmann scalars are given by 



[l-z 



(1 - z^V" 

I ^ i?"^^^i?„^Aa = 48g4^ ^r^. (3.5) 

From these expressions we can see that the spacetime is singular on the hypersurface t = 0. On the other hand, 
although the metric is singular on z = 1, the spacetime is not. Thus, to have a geodesically complete spacetime, 
we need to extend the metric beyond this surface. In order to do so, it is found convenient to study the two cases 
< 2q'^ < 1 and 2q'^ > 1 separately. 



A. < 2?^ < 1 

In this case, introducing two null coordinates u and v via the relations 

r = {-ur - i-vr = f_{u,v), (3.6) 

we find that in terms of u and v the metric and massless scalar field take the form 

ds^ = n''4}'^fl^''-^^'^dudv - fldw^ - f_de\ 

</. = 2g In [/+(«,«)], (3.7) 

where 

n.^>l. (3.8) 

From Eq.(3.6) we can sec that the region t < 0, r > 0, z < 1 in the (t, r)-planc is mapped into the region 
u, V < 0, V > u, which will be referred to as Region // [cf. Fig. 1]. The half line z = 1, f < is mapped to 
V = 0, u <0. The region v > 0, u < 0, which will be referred to as Region 7, is an extended region. Depending on 
the values of n, the nature of the extension is different. In particular, it is analytical only for the case where n is an 
integer. Otherwise, the extension is not analytical, and in some cases the metric and the scalar field even become not 
real in this extended region, as we can see from Eqs.(3.6) and (3.7). To have the extension unique, in the following 
we shall consider only analytical extensions, that is, the cases where n is an integer. Then, from Eqs.(3.7) and (B.7) 
we find that 

(^,„ = -2nq 7 , = -2nq 7 , 

/+ /+ 

i? = ^„,7,.«=nV4i+V"^"'')""' 



2-1/n ( n-i y {uvf 



_ 4i-V"(-^)^"-i _ 4^-V"(-u)^»-i 

„_/-_y2(2n-l)/n ' ^" ^^_^2(2„-l)/„ ' \-^-^> 

Prom these expressions we can see that the spacetime is regular on the symmetry axis v = u'm. Region //, and (j),a is 
always timelike. On the hypersurface v = 0, we have 

<^,„(u,0) = 0, (3.10) 
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and the only non- vanishing component of the energy-momentum tensor T^^, is given by 



(3.11) 



which represents an energy flow, moving from Region II into Region I along the null hypersurfaces u = Const. The 
expansion, Qi, of the null geodesies along the hypersurfaces u = Const, is always positive in this region, while the 
expansion, 0„, of the null geodesies along the hypersurfaces v = Const, is always negative. However, G; becomes 
zero on the hypersurface v = and then negative in the extended region, 7, where v > 0, while G„ is negative even 
in this extended region. Thus, all the cylinders of constant t and r are trapped in the extended region, but not in 
Region //. Then, the hypersurface v = defines an apparent horizon [18,19]. 

It should be noted that the above analysis is very important when we consider boundary conditions on the apparent 
horizon in the next section, as it shows clearly that it is the component (/> „ that represents the energy flow of the 
scalar waves that moves from Region II into Region I, while the component represents the energy flow of the 
scalar field that moves in the opposite direction. Since now Region 7 is a trapped region and no radiation is able to 
escape from this region. This can be seen clearly from Eq.(3.10). 

The singularity behavior in Region 7 depends on the values of n. In particular, when n is an odd integer, from 
Eq.(3.9) we can sec that the spacctime becomes singular on the hypersurface r = or m = —v, which services as the 
up boundary of the spacetime, and the corresponding Penrose diagram is that of Fig. 1. Thus, in this case Region I 
can be considered as the interior of a black hole, and the corresponding solutions represent gravitational collapse of a 
massless scalar field in Region 77. The collapse always forms a black hole. Note that in this case (j)^a is continuously 
timelike in the trapped region, 7, as we can see from Eq.(3.9). 



FIG. 1. The Penrose diagram for the solutions given by Eqs. (3.6) and (3.7) with n > 2 being an integer. The cyhnders of 
constant t and r are all trapped in Region I where Qi < and @i@n > 0, but not in Region //, where 0; > and QiQn < 0. 
The dashed line v = represents the apparent horizon. When n is an odd integer, the spacetime is singular on the horizontal 
double line r = 0, and when n is an even integer, the spacetime has no curvature singularity there. 

When n is an even integer, from Eq.(3.9) we can see that the Ricci and Kretschmann scalars are all finite at r = 
in Region 7, but both of 6/ and B„ become singular there. Thus, anything that moves along the null geodesies, 
defined by or n^, will be crashed to zero volume by the infinitely large contraction. Then, the hypersurface r = 
now represents a topological boundary of the spacetime, and the corresponding Penrose diagram is also given by Fig. 
1, but now the spacetime is free of curvature singularities on the double horizontal line r = 0. 

It should be noted that apparent horizons and black holes are usually defined in asymptotically flat spacetimes 
[18]. To be distinguishable, Hayward called such apparent horizons as trapping horizons and defined black holes by 
the future outer trapping horizons [19]. For the sake of simplicity and without causing any confusions, in this paper 
we shall continuously use the notions of apparent horizons in the places of Hayward's trapping horizons, and define 
black holes in a little bit more general sense than that of Hayward in non- asymptotically flat spacetimes. 



r = 



o 




B. > 1 



In this case, introducing the two null coordinates u and v via the relations 
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t = u + V, r = V — u, 



(3.12) 



we find that the metric and scalar field are given by 

29^ 



{u + vY 



uv 

2q\n[-{u + v)]. 



dudv — (u + v^dw^ — {u — v^dO"^, 



(3.13) 



To study the physics of the spacetime near the hypersurface v = 
null geodesies along the hypersurface u = Const., say, u = uq, 

V — 2q — -V 

V[Uo + V) 



or 2 = 1 in some details, let us consider the radial 



0, 



(3.14) 



where an over-dot denotes the ordinary differentiation with respect to the affine parameter A along the null geodesies. 
Then, near the hypersurface v = 0, Eq.(3.14) has the solution 



(&1A + 62) 

>iA+62 



-i/(2g" 



2q2 



> 1, 

= 1, 



(3.15) 



■AS V ^ 0, we must have A ±00. That is, the "distance" 
say, (u.v) = {uo,vq < 0), along the null geodesies u = uq is 



where bi and 62 arc the integration constants. Thus, 
between the point (mo,0) and any of the other points, 

infinite. Therefore, when 2q'^ > 1 the hypersurface v = actually represents a natural boundary of the spacetime, 
and there is no need to extend the solutions beyond this surface, since now Region II is already geodesically maximal. 
It should be noted that, although there is no spacetime singularity on the half-line v = O.u < 0, the spacetime is 
singular at the point {u,v) = (0,0), as can be seen from Eq.(3.9). The corresponding Penrose diagram is given by 
Fig. 2. 



(0, 0) 




FIG. 2. The Penrose diagram for the solutions given by Eq. (3.13) for 2q^ > 1. The spacetime is geodesically maximal in 
the whole region u,v < 0, v > u, and singular only at {u,v) = (0,0). The two surfaces of constant t and r are not trapped, 
because now we always have 0( > and Qn < 0. The only exception is on the surface v = 0, u < where Qi{u, 0) = and 
6„(w, 0) < 0. The dashed lines represent the hypersurfaces <f){u,v) = Const., which are always spacelike. 
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On the other hand, from Eq.(B.7) we find that 



uv 



{u + vf 



UV 



2q^ 



(3.16) 



Thus, in the whole spacetime now we always have 9; > and 0„ < 0, except on the half-hypcrsurface v = 0, u < 
where we have Qi = 0, 6„ < 0. That is, all the two surfaces of constant t and r are not trapped for u,v < 0, and 
become marginally trapped only on the half surface v = 0,u < 0. In addition, we have 



R = 



(3.17) 



which is always positive for u,v < 0, and zero only when v = 0. That is, the scalar field is always timelike, except 
on the hypersurface v = where it becomes null. Thus, in this case the corresponding solution can be considered as 
representing gravitational collapse of a massless scalar field. Although now no black holes are formed, a point-like 
spacetime singularity is indeed developed at the point («, v) = (0, 0), as we can see from Eq.(3.5). It is interesting to 
note that this singularity is not naked, and an observer can see it only when he/she arrives at that point. 



IV. LINEAR PERTURBATIONS OF THE SELF-SIMILAR SOLUTIONS 



To sec if the above solutions represent critical collapse, we need to do their linear perturbations, because by 
definition a critical solution has one and only one unstable mode. To study such perturbations, it is found convenient 
to use the self-similar variables r and z defined by Eq.(A.l) but still work in the (f, r)-coordinates. Then, the linear 
perturbations can be written as 



F{T,z) = Fo{z) + eFi{z)e 



kT 



(4.1) 



where F = {M, S, V, (/?}, and e is a very small real constant. Quantities with subscripts "1" denote perturbations, and 
those with "0" denote the background self-similar solutions given by Eq.(3.4). It is understood that there may be 
many perturbation modes for different values (possibly complex) of the constant k. Then, the general perturbations 
will be the sum of these individual ones. Modes with Re{k) > grow as r — > oo and are referred to as unstable 
modes, and the ones with Re{k) < decay and are referred to as stable modes. 

It should be noted that in writing Eq.(4.1), we have already used some of the gauge freedom to write the perturba- 
tions such that they preserve the form of the metric (2.1). However, this does not completely fix the gauge. We shall 
return to this point later when we consider the gauge modes. 

To the first order of e, the Ricci tensor is given by Eqs.(A.4)-(A.10). Applying them to the background solutions 
given by Eq.(3.4), and using the Einstein field equations (3.1), we find that there are only four independent equations, 
which can be cast in the form, 



kMi{z) = 2z^S'l + 2zVl + 2 



+fc 1 - 4q2 



1-z 



2 I Si+ ^qzip'-y. 



S[ + kVi 



(4.2) 



and 



z{l- z'^) VC + 2k)z'^] Vl - k^zVi = kzSi - (l - z'^) S[, (4.3) 

z{l- z^) ip'l + [1 - (1 + 2k)z^] ip[ - k^zifii = 2qz {zS[ + kSi) , (4.4) 

z{l- z"^) S'l + 2 {I- kz'^) S[ + k{l - k)zSi = 0, (4.5) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. It can be shown that 
Eq.(4.5) has the general solution, 

5i {z) = 1 [ci (1 + ^)2-'= + C2 (1 - zf-'^] , (4.6) 
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where ci and C2 are two integration constants. Substituting the above solution into Eqs.(4.3) and (4.4) we find that 
these two equations can be written in the form 



J2 7. J7. 1 



fi{y), 



(4.7) 



where y = z^, {Zi} = {Vi,ipi}, a = b = k/2, e = 1, and 
h{y) = kzSi-{l-z^) S[ 



jci [1 - (2 - fc)yi/2 + y] (1 + + c, [l + {2 - k)y'/^ + y] (l - y'^'Y''} 



f2{z) = 2qz{zS[+kSi) 



= -2q jci [(1 - fc) - + ' V C2 [(1 - fc) + y'/^] (l 



.1/2 



1-fe 



(4.8) 



Eq.(4.7) is the inhomogeneous hypergeometric equation [20], and the general solution of the associated homogeneous 
equation is a linear combination of the two independent solutions, ^^''(z), where 



Fi'\z) = F(^k,^k;l;z' 



1. 1 



Ff)(z) 



Flh,^k;k;l 



(4.9) 



with F{a, b; e; z) denoting the hypergeometric function. From the above two independent solutions, we can construct 
particular solutions of the inhomogeneous equation (4.7), and then find that the general solutions for Vi{z) and fiiz) 
can be written as 



V^iz) = (af ) + Af\z)) )(.) + (aW - A«(.)) )(.), 

^i(^) = (4^^ + (^)) Fi'\z) + (a« - aW(.)) )(.), 



(4.10) 



where a^'^'s are integration constants, and 



fAz)Fl'\z)dz 
z{l- z'^)A{z)' 



Aiz) ^ Fi'\z)l (f«(.)) - Fi'\z)l (Ff )(.)) . 



(4.11) 



To have physically acceptable perturbations, we need to impose boundary conditions. In General Relativity, this 
is a very subtle problem and there are no fixed rules to follow. In this paper we shall choose the axis r = and the 
hypersurface 2; = 1 as the places where we impose the boundary conditions. 

Since the axis for the background solutions is regular, and the conditions (2.2)-(2.3) and (2.9) are satisfied by them, 
we would expect that the linear perturbations also satisfy these conditions. In particular, it can be shown that the 
condition (2.9) requires Mi(0) = 0, and the ones (2.2) and (2.3) require, respectively, Gi{z) and zG'i{z) 0, as 
z ^ 0, where Gi{z) = Si{z) + Vi{z). On the other hand, the condition that R = 4>,a4''" is regular on the axis further 
requires zip'i + k(pi — > ^. In summary, on the axis we shall impose the following conditions, 

Mi(z) ^0, 
Gi(^) -^0, 
zG'i(z)^0, 

z^\{z) + kipi{z) ^ 0, (4.12) 



^It should bo noted that this condition is not independent of the ones (2.2) and (2.3). In fact, using the Einstein field 
equations (3.1) we can deduce it from Eqs.(2.2) and (2.3). However, without loss of generality, in this paper we shall impose it 
independently. 
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as z ^ 0. 

On the other hand, when < 2<f' < 1 the hypersurface z = 1 is an apparent horizon [cf. Fig.l ], and we required 
that the background solutions be analytical across it with respect to the null coordinate v. Otherwise, it was found 
that the extension was not unique. Clearly, this condition should hold also for the perturbations. In addition, since 
the hypersurface z = 1 represent an apparent horizon and Region J is a trapped region, so nothing should be able 
to escape from it. In particular, for the scalar field this implies (j>^v{u,0) = [cf. Eq.(3.11)]. On the other hand, as 
shown in [21], the gravitational wave component that moves out of Region / is represented by \E'o, which is a function 
of Vvv, Vv, M^y and S^y. Then, we can see that the condition that no gravitational waves come out from Region I 
requires V,vv, V,v, S,v, M^,., — > as u — > 0. Changing to the self-similar variable z, it can be shown that these 
conditions are equivalent to 

Fi{z) ^ analytical with respect to v, 

(l_^)(n-l)/n^^^O, (4.13) 

dz 

(l_^)(n-2)/«^}^^0^ (0<2g2<l), 

as z ^ 1, where Fi{z) = {Mi, Si,Vi,ipi}. 

When 2q^ > 1 the hypersurface z = 1 is a future null infinity of the spacetime and marginally trapped [cf. Fig. 
2]. Since it is the future null infinity, we would expect that the perturbations be finite there. On the other hand, 
the hypersurface z = 1 is also marginally trapped, so there should have only outgoing scalar field and gravitational 
waves. Note that the amplitude of gravitational wave components is always proportional to ^ (1 — z)'^'^ [21]. 
Thus, now we must require (1 — z)^"? (Vi_vv, Vi.v, Si,y, Mi^^) — > 0, as z — > 1. On the other hand, using the Einstein 
field equations we can show that these conditions also imply (1 — z)^'^(^i,j, 0. In terms of the self-similar variable 
z, and noticing that now t and r are given by Eq.(3.12), we find that these conditions can be written as 



Fi (z) —^ finite 

Fi( 

dz 

rfz2 



dz 

(l--)^'^^^-O, (2g2>l), (4.14) 



as z — »■ 1. 

Once we have the boundary conditions, let us first consider the gauge modes. We note that the metric (2.1) is 
invariant under the coordinate transformations, 

t = a{t + f)+ b{t- f), 

r = a{t + r) -b{t-f), (4.15) 

where aif + f) and b(t — f) are arbitrary functions of their indicated arguments, subject to a'b' ^ 0. Thus, let us 
consider the gauge transformations 

t^t + e [A{t + r) + B{t - r)] , 

r -*r + e[A{t + r) - B{t -r)], (4.16) 

where A{t + r) and B(t — r) are other arbitrary functions. Then, we find that under the above coordinate transfor- 
mations the resultant perturbations are given by 



with 



F{t,z) = Fo{z) + 6F{t,z), (4.17) 



SM{t, z) = 2{A' + B') - ^^i^ [{A -B) + {A + B)z] , 

5SiT, z) = i [{A -B) + iA + B)z] , 
SV{t, z) = -^[iA-B) + {A + B)z] , 

6ip{T,z) = ^{A + B). (4.18) 
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In order to have the above expressions be in the form of Eq.(4.1), we must choose 



B(t-r)=Ciio \l + zy-'^, (4.19) 
for which Eq.(4.18) can be written as 

6F{T,z)^Fiiz)e''\ (4.20) 

with 

Miiz) = 2 {ci [(1 ~k) + 2q^z] (1 + z)-'' - C2 [(1 - fc) - 2q^z] (1 - z)-'^} , 

S,{z) = - [ci(l + zf-'^ + C2(l - zf-"] , 
z 

V,{z) = [ci(l + zy^ + C2(l - z)-''] , 

z 

M^) = -25 [ci(l + z)'-'' - C2(l - z)'-''] , (4.21) 

where Ci and C2 are arbitrary constants. It can be shown that the above perturbations don't satisfy the boundary 

conditions (4.12)-(4.14), or in other words, the boundary conditions imposed in this paper limit all the gauge modes. 

Now let us consider the above boundary conditions for the particular background given by Eq.(3.4). We first 
consider the conditions at the axis given by Eq.(4.12). Let us first note that 

Fi'\z)^l + ^k'z' + 0{z^), 

f[^\z) = A2{k) |[Ai(fc) -2ln{z)]Fi'\z) + ^kik - 2)z^^ , (4.22) 

as z ^ 0, where 

, ^ 2(2 -k) , T(k) 

and r(fc) denotes the gamma function. Inserting Eqs.(4.22) into Eqs.(4.11) and (4.10), after tedious calculations we 
find that 

zip\{z) + kipi{z) ^ -2fc4^^^2(A;) ln(.2) + Aoik), (4.24) 
as — > 0, where Ao{k) is a finite constant. Thus, the last condition of Eq.(4.12) requires 

= 0. (4.25) 

Similarly, one can show that 

Gi(z) = Siiz) + Viiz) = ^^'''^'''^ - 2a['^A2ik) \n{z) 

z 

+ (2 - k) (ci - C2) + (^a^^'> Ai{k)A2{k) + ^) + O (z) , (4.26) 

from which we can see that the second condition of Eq.(4.12) requires 

ci = -C2 = c, a^^^ = 0, ^ = 2{k - 2)c. (4.27) 
Once Eq.(4.27) holds, it can be shown that 

zG\{z)^0{z^), 

M,{z) ^ O (z'') , (4.28) 

as z ^ 0. That is, the first and third conditions of Eq.(4.12) do not impose further restrictions on the free parameters. 

Now let us turn to consider the boundary conditions at z = 1. It is found convenient to consider the cases 
< 2g'^ < 1 and 2g'^ > 1 separately. 
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A. < < 1 



In this case we have 1 — z as ?; ^ 0. Then, we find that 

Si{z) ~ c (22-'= - (-v)(2-fc)»j , (4.29) 
as V — > 0. Thus, in order to have be analytical, the constant k has to take the values, 

k = 2 , (to>1), (4.30) 



n 



where m is a positive integer (m = 1, 2, 3, ...). Since when A; < the corresponding modes are stable, which we are not 
so interested in. Thus, in the following we shall consider only the case where fc > 0, which together with Eq.(4.30) 
implies 2n > m > 1. To study the boundary conditions &t z = 1 further, let us consider the cases k ^ 1 and k = 1 
separately. 

Case A) 7^ 1: In this case it can be shown that 

f[^\z) = B,{k)F^l\z) + B^{k)F[^{z)x^-\ 

Fi'\z) = Fil\z), (Ml), (4.31) 



where 



fII\z) = l + ^kx + Di{k)x'^ + O (x^) , 

F^^\z) = 1 + h2 - k)x + D2 {k)x^ + i?3 (A;)a;3 + O (x^) , (4.32) 



with x = l — z'^ and 



^i(^)-r2 (i-ifc)' ^^(^)-T^' 

Inserting the above expressions into Eqs.(4.2), (4.10) and (4.11) and after tedious calculations we find that 

kMi{z) = ifc {c(fc - 2) (fc - 1 - 2q^) 2^ - 2Bj,{k) {k - l)af ^ + 2fcg4^^] } x^'*" 
-2 {c(fc - 2) (fc - 1 - 2g2) 2'= + 2(1 - k)B2{k) af ^ + 2g4^'] } x''' 

+^3(0, k,q)+0 (x'=) + O (,t2-'=) , (4.34) 
(^i(z) = [S2(fc).xi-'= + Bi{k)] +0{x), (4.35) 
Vi{z) = a'i^ B2{k)x^-'' + O (x^-'^) , (fc 1), (4.36) 

(2) 

as z ^ 1, where ^3(0, fc, q) is a finite constant, and a\ is given by Eq.(4.27). From the last two equations we can see 
that the analytical conditions of Vi{z) and ipi{z) require fc < 1, which together with Eq.(4.30) implies 

2n> m>n. (4.37) 

The condition that Mi (z) is analytical across the hypersurface z = 1 requires 



^ 2 V 2(fc - l)B2(fc) 

On the other hand, it can be also shown that 
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„1-1 



/n dM^jx) c{k - 2) 



X 



dx 
dVi{x) 

dx 

i_i/„d^i(^ 



[(fc - 1 - 2q^) 2^ + 4(1 - k)B2{k)\ a;(™-"-i)/" + O (a;^"-^)/") , 



l-l/n 



(1 - A:)al')B2(fc)x('"-"-2)/" + O (x^-^/") + O (x^^-^)/") , 



as a; — » 0. Thus, the last two conditions of Eq.(4.13) further require 

2n>m>n + 2, {k 1). 
Case B) = 1: In this case it can be shown that 

F^'\z) ^ F^l\x) - Fil\x)ln{x), 



(4.39) 



(4.40) 



but now with 





_ 2 

TT 




_ 1 

TT 

oo 


Co 


n=0 



Ff)(z)=^F«(x), (fc=l), 



Co + J(Co - l)x + Y^(6Co - 7)a;2 + O (x^) 



(4.41) 



1 + ia;- 
4 



64 



0(.3) 



(n + l)(n + 2)' 



Then, it can be shown that 



Vi{z) 



2c 



ln(a;) +0(.Tln(a;)) , (fe = 1), 



(4.42) 



(4.43) 



as a; ^ 0. Thus, the analytical condition of Vi{z) across the hypersurface a; = (2 = 1) requires c = 0. It can be 
further shown that 



>pi{z) ln(a;) + O (xln(a;)) , {k = 1), 



(4.44) 



(2) 

as X ^ 0. Thus, the analytical condition of (pi{z) across the hypersurface x = requires =0. Considering 
Eqs.(4.25) and (4.27) we find that the boundary conditions in the present case limit all the perturbations, that is, 

M,{z) = Sr{z) = V,{z) = Mz)=0, {k = l). (4.45) 
Therefore, from Eq.(4.40) we can see that for any given n, the solution has 

N = n-3, (0<2g2<i), (4.46) 

unstable modes. In particular, the solution with n = 2 or n = 3 has no unstable mode, and consequently is stable 
with respect to the linear perturbations. The solution with n = 4 has only one imstablc mode, which may represent 
a critical solution, sitting on a boundary that separates two attractive basins in the phase space. All the solutions 
with n > 4 have more than one unstable modes and are not stable with respect to the linear perturbations. 



B. 2q^ > 1 

In this case the boundary conditions at 2; = 1 are those given by Eq.(4.14). From Eq.(4.6) we can see that to have 
Si{z) be finite as 2; ^ 1, we must assume that Re{z) < 2. To study the boundary conditions further, let us first 
consider the case k ^ 1. 
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When k ^ 1, it can be shown that Eqs.(4.31)-(4.36) also hold in the present case. Thus, the conditions that Vi{z) 
and tpi (z) are finite as 2; ^ 1 require that 

Re{k) < 1, (4.47) 

while the condition that Mi{z) is finite further requires that the constant a2 has to take the values given by Eq.(4.38). 
On the other hand, from Eqs.(4.34)-(4.36) we can see that the last two conditions of Eq.(4.14) require 

Re{k) < 2(t - 1. (4.48) 

When k — \,\\, can be shown that Eqs.(4.41)-(4.44) hold in this case, too, and following the analysis given there 
we find that the boundary conditions at z = 1 limit all the perturbations for fc = 1. 

Therefore, for any given q with > 1, Eqs.(4.47) and (4.48) show that there always exists a continuous spectrum 
of k such that 

0<i?e(fc)<i?e(A;™„.), (2^^ > 1), (4.49) 

where 

i?e(fc_.) = {2,2_i, f<,2'<i. (4.50) 

That is, in the case > 1, there are infinite numbers of unstable modes. However, when = 1, from Eq.(4.48) 
we find that 

Re{k) < 0, (2g2 = 1). (4.51) 
Thus, the solution with 2q'^ = 1 is stable against the linear perturbations. 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper we have first introduced the notion of homothetic self-similarity to four-dimensional spacetimes with 
cylindrical symmetry, and then presented a class of exact solutions to the Einstein-massless scalar field equations, 
which is parameterized by a constant, q. It has been shown that for < 2g^ < 1, the corresponding spacetimes have 
black hole structures but with cylindrical symmetry. These black holes are formed from the gravitational collapse of 
a massless scalar field. When 2q^ > 1 the corresponding solutions also represent gravitational collapse of the scalar 
field but no black holes are formed. Instead, a point-like singularity is developed, which is not naked and can be seen 
by an observer only when he/she arrives at the singularity. 

Then, the linear perturbations of all these solutions have been given analytically in closed form in terms of hy- 
pergeometric functions. After properly imposing boundary conditions at the axis and on the horizons, it has been 
shown that the solutions with n = 2, 3 and the one with 2q'^ = 1 are stable, where n is an integer and given by 
n = 1/(1 — 2(7^). For any given n > 4, the corresponding solution has TV = n — 3 unstable modes. In particular, 
the one with n = 4 has precisely one unstable mode, which may represent a critical solution sitting on a boundary 
that separates two attractive basins in the phase space. The solution for any given q with 2q^ > has a continuous 
spectrum of unstable eigen- modes, given by Eqs.(4.49) and (4.50). 

It should be noted that in this paper we have shown that black holes can be formed from gravitational collapse of the 
massless scalar field, and the ones with n = 2, 3 are stable against the linear perturbations. However, these spacetimes 
are not asymptotically flat in the radial direction, and thus may not be considered as representing counter-examples 
to the hoop conjecture [22]. To have an asymptotically flat spacetime in the radial direction, we may restrict the 
distribution of the scalar field only to a finite region, say, r < ro{t), and then join it with an asymptotically flat region. 
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APPENDIX A: THE RICCI TENSOR AND ITS LINEAR PERTURBATIONS IN TERMS OF 

SELF-SIMILAR VARIABLES 

The general metric for cylindrical spacctimes with two hypcrsurfacc orthogonal Killing vectors takes the form of 
Eq.(2.1). Introducing the self-similar dimensionless variables r and z vis the relations, 

where is a dimensional constant, we find that the non- vanishing components of the Ricci tensor are given by, 

= ^ {2^^5,.. - (1 - z") (zM,,, + - (^,,2 ^ y 2^ + z (1 + z") M^,S,, + \z'S^, 

+Z [25 + M^rr - {S,r^ + Vr^) +{2 + M^) S .r + M,^] 

+z^ [45,, + 2M^r. - 2 (5,5, + VrV,) + (M,,5, + M,,5,)]} , 
Rtr = ^ {2^'5 - z^ (5,2 + V,^) -2z{l- ^5 ,) M,, + 4^5 , 

+z [25,, - (5,5 , + + (M,,5 , + M,,5,)] + 2 (5, - M,,)} , 



Rrr = ^ {2z5 „ + z (1 - 2^) M,„ - z (5,,' + V/) +z{l + z^) M,,S,, - 2 (l + z^) m,, + 45,, 

-z [M,„ + (1 - 5,) M,,] - [2M,,, - (M,,5, + M,,5,)]} , 
i?22 = ie*^+^-^ (1 - [(5„ - y,,) - 5 , (5, - F,)] + 2z (1 - _ -I/,) - 2 (1 - z5,,) 

-z3 [2 (5,, - V,,) - 25,5, + (5,y, + 5,y,)] 
-z^ [(5„ - y„) + (1 - 5,) (5, - y,)]} , 

i?33 = ie^-^-^ {z^ (1 - z2) [(5 „ + - 5, (5 , + y,)] + 2z (1 - z^) (5^ + y^) _ 2 (1 - z5,) 

-z^ [2 (5„ + - 25,5, - (5,V, + 5,y,)] 

-z2 [(5 „ + V,,) + (1 - 5 ,) (5 , + V,)] } . (A.2) 

On the other hand, it can be shown that the Klein-Gordon equation, □</> = 0, for the massless scalar field takes the 
form 

Z(l - z'^)4>^zz - 2z^4>^rz - Z(t),rr + [z{zS,;, - 1) + zS^r] 'P,r 

+ [(1 - z^){2 - z5,,) + z^S,r] ct>,z = 0. (A.3) 

Now let us consider the linear perturbations of Eq.(4.1). To first order in e, it can be shown that the non- vanishing 
components of the Ricci tensor are given by 

ijW = [2z^s'; - z (1 - z2) M'{ - 2zVoV/ - [2 - 2(1 + k)z'' - z (l + z^) S'^] M[ 

+z [4(1 + k)z +{1 + z^) M'o - 2z^S'q] S[ - 2kz^V^Vi +kz{l + k + zS'^) Mi 

+fcz[2(l + A;) + z(M^-25^)]5i}, (A.4) 

(2+A;)t 



-2z (1 - zS'q) M[ - kzV^Vi -k{2- zS'o) Mi + k[2 + z [M'^ - S'q)] Si) , (A.5) 

p(2+fe)-!- 

= —-2- {2z5;' + z (1 - z2) M{' - 2zV^Vi + [4 + z (l + z^) - 2z5^] S[ 

- [(1 + z2) (2 - z5^) + 2kz'^] M[-kz{l + k- zS'q) Mi + kz^M^^Si} , (A.6) 

4^2^ = e'^^ I (Ml + Fi - 5i) 4? + ize^''+^''-^°i?(;,) 1 , (A.7) 



i?« = e'^- \ (Ml - Fi - 5i) + ize^-^o-so^(i) ^ ^ (A.8) 
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where R^^ and r!^^ are the corresponding components of the Ricci tensor for the background solution, and 

=z{l- {S'l - Vi') + [2kz^ + (1 - z^) {zS'o - 2)] VI + [4 - 2(1 + k)z'' - z {l - z^) {2S'o - V^)] S[ 
+kz {1 + k- zS'q) Vi - kz [(1 + k) - z (25^ - V^)] Si, (A.9) 

=z{l- z^) {S'l + O - [2kz^ + (1 - z^) {zS'„ - 2)] VI + [4 - 2(1 + k)z^ - z {l - z^) (25^ + V^)] S[ 
-kz {1 + k- zS'o) Vi - kz [(1 + k)-z {2S'o + ^o')] ^i- (A.IO) 



APPENDIX B: APPARENT HORIZONS IN SPACETIMES WITH CYLINDRICAL SYMMETRY 

In [23], the ingoing and outgoing radial null geodesies were studied in double null coordinates, and the corresponding 
expansions of them were calculated. In this paper, we shall define apparent horizons in spacetimes with cylindrical 

symmetry, using those quantities. Before doing so, we would like first to note the difference between the double null 
coordinates used in this paper and the ones used in [23]. As a matter of fact, the roles of u and v are exchanged in 
this paper. 

To write the metric (2.1) in the double null coordinates, let us first introduce the two null coordinates u and v via 
the relations, 

t = a(u) +/3(u), r = a{u) - j3{v), (B.l) 
where a{u) and (}(v) are two arbitrary functions of their indicated arguments, subject to 

a'{u)l3'{v) + 0, (B.2) 
where a prime denotes the ordinary differentiation. Then, in terms of u and the metric (2.1) takes the form 

ds^ = 2e^^^"'^^dudv - r^e-^^*''-) (e^^*''-)^^;^ + e-^^''''^d9^) , (B.3) 

where 

a{u, v) = ^ [In {2a' 13') - M] . (B.4) 



Introducing two null vectors l'^ and n"^ by 
we find that 



«A = ^='^A, n,^—^=5l (B.5) 



= = n^;.n", (B.6) 

which means that each of them defines an afhnely parametrized null geodesic congruence. In particular, defines the 
one moving along the null hypcrsurfaccs u = Const., while defines the one moving along the null hypersurfaces 
V = Const. Then, the expansions of these null geodesies are defined as, 

e„ = <?"V;/3 = e-2-^, (B.7) 

where 

Tl=\\d^-d^\\-\\de-de\\ = r^e-^. (B.8) 

We call the cylinders of constant t and r trapped if 0/0,, > 0, marginally trapped if 0;0,i = 0, and untrapped if 
0(0„ < 0. An apparent horizon is defined as a hypersurface foliated by marginally trapped surfaces [18,19]. 
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